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INTRODUCTION 


Cross  derivatives  are  generated  in  the  vorticity-transport  equation 
when  it  is  transformed  to  non-orthogonal  coordinates.  The  presence  of 
cross  derivatives  in  the  equation  cause  some  difficulties  in  the  method 
of  Ref.  [1].  More  specifically,  they  arise  in  the  integration  of  the 
vorticity  diffusion  form  over  the  control  volumes.  These  difficulties 
are  eliminated  by  modifying  the  diffusion  form  of  the  equation  and 
the  Integrating  factor  K.  This  treatment  is  suitable  for  both 
elliptic  and  parabolic  type  partial  differential  equations.  Details 
of  the  derivation  are  given  in  the  following  sections. 
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ANALYSIS 

The  vorticity-Cransport  equation  in  non-orthogonal  coordinates  can  be 
written  in  the  following  form: 


aoi  +  bio  +  cui  -  du>  -  eto  *  0  ,  (1) 

xx  xy  yy  x  y 

where  <o  is  the  vorticity,  a,  b,  c,  d  and  e  are  functions  of  x  and  y.  For  both 

-2 

elliptic  and  parabolic  partial  differential  equations,  b  -  4ac  <  0. 

Direct  application  of  the  method  of  Ref.  [1]  to  Eq.  (1)  yields  the 
following  diffusion  form: 

(Kaui  )  +  (Kcm  )  +  Kbm  -  0  ,  (2) 

x  y  y  y 

where 

a  e 

K  * - x - y  ,  a  and  c  *  0  .  (2a) 

a  c 

a,  b,  c,  d  and  e  are  the  linearized  values  of  the  coefficients  a,  b,  c,  d 
and  e  respectively. 

The  Integration  of  the  third  term  in  Eq.  (2)  over  the  control  volume  of 
Fig.  1  does  not  yield  a  suitable  form  for  the  5-point  and  9-point  schemes  [1]. 

The  diffusion  form,  Eq.  (2),  is  modified  to  eliminate  the  integration 
difficulties  encountered  in  the  cross  derivatives  term. 

The  first  and  second  order  derivatives  with  respect  to  x  are  combined 
together  in  one  term  using  the  integrating  factor  in  the  x-directlon. 

E*P  (7  *)  [E*P  (“  f  *)  +  b“xy  +  c"yy  -  e«y  -  0  .  (3a) 


2 
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Equation  (3a)  Is  rewritten  as 


EXP  X)  [EXP  ("  J  X)  (a“x  +  b“y)]  +  7  b“y  +  C"yy  "  e“y  "  0  *  (3b> 

The  last  three  terms  In  Eq.  (3b)  are  combined  together  In  one  term  using 
the  Integrating  factor  in  the  y-dlrection. 


EXP  (j  X)  [EXP  ("  7  X)  (»“x  +  b“y)]  + 

7  X 


Exp  (--a-^--b-  y)  [Exp  (-  y)  CM  ]  -  0  . 

y 

By  multiplying  and  dividing  the  first  term  in  Eq.  (3c)  by 
Exp  [-  — - —  y)  and  the  second  term  by  Exp  (-  —  x),  one  obtains 


[K(ao>x  +  b»  )]  +  [Kcw  ]  -  0  , 

7  x  7  y 


where 


„  d  ea  -  db 

K  - - x - y 

a  ac 


The  modified  diffusion  form,  Eq.  (4),  is  then  integrated  over  the  control 
volume  shown  in  Fig.  1  to  obtain  Eq.  (5). 

y  =  Ay/2  y  -  Ay/2 

/  K(au)x  +  bajy)j  dy  -  /  K(au>x  +  bwy)J  dy 


y  ”  -  Ay/2 


x  ■  Ax/2 


x  -  Ax/2  y  -  -  Ay/2 


x  «  Ax/2 


x  -  -  Ax/2 


+  /  Kcu>y|  dx  -  /  Kc<uy|  dx  »  0  .  (5) 

x  -  -  Ax/2  y  -  Ay/2  x  -  -  Ax/2  y  -  -  Ay/2 
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The  Modified  5-Polnt  Scheme 

The  vortlclty  derivatives  &i  and  ui  at  x  *  ±  Ax/2  and  y  *  ±  Ay/ 2  are 

x  y 

approximated  using  central  finite  difference  equations.  The  vorticity 
derivatives  at  x  -  ±  Ax/2  are  approximated  as  follows: 


(“i+l.j+1  +  “i.j+1  “  “i+l.J-1  "  “i.j-JMAy 


-  Ax/ 2 


+  “i-i.j+1  '  “i.1-1  •  Vi,j-I>4iy  •  (6b) 


Upon  substitution  into  the  integrated  form  of  Eq.  (4),  i.e.,  Eq.  (5), 


one  obtains 


Mi+u  +  Mi-u  +  +  Mi,j-i +  (B“W,j+i 


+  Mi+i.j-i  +  (B“)i-I,j+1  +  Mj-u-i  -  (B“)i,j  (7) 

where 

B1+1^  *  a  [Exp  (-  p  -  q)  -  Exp  (-  p  +  q)/(2Ax2q)  (7a) 

Bi,.1+1  *  C  ^EXP  ("  p  _  <0  ~  ExP  (p  ~  9)]/(2Ay2p) 

+  b  [Exp  (-  p)  -  Exp  (p)]  [Exp  (-  q)  -  Exp  (q)/(8AxAyq)  (7b) 

Bi+1 ,1+1  *  b  fExp  (“  P  -  q)  -  Exp  (-  p  +  q)]/[8AxAyq)  (7c) 

j  -  a  [Exp  (-  p)  +  Exp  (p)]  [Exp  (-  q)  -  Exp  (q)]/(2Ax2q) 

+  c  [Exp  (-  q)  +  Exp  (q)  ]  [Exp  (-  p)  -  Exp  (p)]/(2Ay2p)  (7d) 
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dAx 

p  "  TT 


(7e) 


q 


(ea  -  db)  Ay 
2ac 


(7f) 


The  rest  of  the  coefficients  in  the  above  Eq.  (7)  can  be  obtained  by 
permutations  of  (7a)  through  (7 f ) . 


The  Modified  9-Point  Scheme 

The  vorticity  flux  components  and  a>y  are  approximated  on  the  control 
volume  surfaces  using  first  order  polynomials  in  x  and  y  [1],  Also,  the 
vorticity  derivatives  oiy  at  x  *  i  Ax/2  are  approximated  using  Eqs.  (6a)  and 
(6b),  Upon  substitution  in  the  integrated  form,  Eq.  (5),  we  obtain  the 
following  equation: 

Mi+i,j  +  +  Mi,j+i  +  Mi,j-i  +  Mi+1,.1+1 

+  (^i+lj-l  +  ^^i-l.j+l  +  ^A“^i-l,j-l  “  ^Aa>^i, j  (8) 

where 

A1+  -  a  Exp  (-  p)  -8^-  +  2c  cosh(q)  (8a) 

’’  Ax  Ay 


i.J+1 


c  Exp  (-  q)  +  2a  cosh(p)  + 

Ay  Ax 


b  [Exp  (-  p)  -  Exp  (p)  ]  [Exp  [-  q)  -  Exp  (q )  ]/ (8AxAyq )  (8b) 
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“  Ai+i  i+i  •  3  Exp  C"  p)  ^  +  c  Exp  ^  ^ 

x+l  ,.1  +  1  AX  Ay 

-  b  [Exp  (-  p  -  q)  -  Exp  (-  P  +  q )  ]/  (8  AxAyq ) 

A  =2  [cosh(q)  +  cosh(p)  ^p-] 

1 » i  Ay  Ax 

f(t)  =  [(l  +  t)  Exp  (-  t)  -  l]/(2t)2 

g(t)  *■  [[1  "  t)  Exp  (-  t)  +  [1  +  t)Exp(t)  -  2]/[2t)2 
g(t)  is  an  even  function  in  t;  g(t)  *  g(~  t)  >  0  and  |g(t)| 


(8c) 

(8d) 

(8e) 

(8f) 

|f[t)|  for 


for  all  values  of  t. 
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U) 
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NOMENCLATURE 

Coefficients  defined  in  Eq .  (1). 

Linearized  coefficients  a,  b,  c,  d  and  e, 
respective ly . 

Functions  defined  in  Eqs.  (8e)  and  (8f). 

Function  defined  in  Eqs.  (2)  and  (4). 

Parameters  defined  in  Eqs.  (7e)  and  (7f). 
Non-orthogonal  coordinates. 

Vorticity . 

Grid  spacing  in  x  and  y  directions,  respectively. 


Subscripts 

x,  y  Refer  to  derivatives  in  x  and  y  directions, 

respective ly . 

i,  .1  Refer  to  grid  lines  in  x  and  y  directions, 

respectively. 
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Figure 


1.  Schematic  diagram  for  grid  geometry 
and  control  volume. 
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